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Abstract 


We present a detailed and systematic analysis of the nonperturbative, 
nonequilibrium dynamics of a quantum field in the reheating phase of in¬ 
flationary cosmology, including full back reactions of the quantum field on 
the curved spacetime, as well as the fluctuations on the mean field. We use 
the 0(A) field theory with unbroken symmetry in a spatially flat Friedmann- 
Robertson-Walker (FRW) universe to study the dynamics of the inflaton in 
the post-inflation, preheating stage. Oscillations of the inflaton’s zero mode 
induce parametric amplification of quantum fluctuations, resulting in a rapid 
transfer of energy to the inhomogeneous modes of the inflaton field. The 
large-amplitude oscillations of the mean field, as well as stimulated emission 
effects require a nonperturbative formulation of the quantum dynamics, while 
the nonequilibrium evolution requires a statistical field theory treatment. We 
adopt the coupled nonperturbative equations for the mean field and vari¬ 
ance derived in a preceding paper |l] by means of a two-particle-irreducible 
(2PI), closed-time-path (CTP) effective action for curved spacetimes while 
specialized to a dynamical FRW background, up to leading order in the 1/A 
expansion. Adiabatic regularization is employed to yield a covariantly con¬ 
served, renormalized energy-momentum tensor. The renormalized dynamical 
equations are evolved numerically from initial data which are generic to the 
end state of slow roll in many inflationary cosmological scenarios. The initial 
conditions consist of a large-amplitude, quasiclassical, oscillating mean field 
(4>) with variance (<b 2 ) — (<£>) 2 around the de Sitter-invariant vacuum. We 
find that for sufficiently large initial mean-field amplitudes > M P /300 (where 
M P is the Planck mass) in this model, the parametric resonance effect alone 
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(in a collisionless approximation) is not an efficient means to “preheat” the 
quantum field. For small initial mean-field amplitude, damping of the mean 
field via parametric amplification of quantum fluctuations is seen to occur, 
and in this case can be adequately described by prior analytic studies with 
approximations based on field theory in Minkowski spacetime. Our results in¬ 
dicate that the self-consistent dynamics of spacetime plays an important role 
in determining the physics of the post-inflationary Universe. This study calls 
into question the validity of general claims made without full consideration of 
the self-consistent dynamics of spacetime and quantum fields. 

PACS number(s): 98.80.Cq, 04.62,+v, 05.70.Ln, 11.15.Pg 
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I. INTRODUCTION AND SUMMARY 


The inflationary Universe P |T3|| has for over a decade been the new paradigm for ad¬ 
dressing many basic issues in cosmology such as the flatness-entropy problem, the horizon- 
homogeneity problem, and the fluctuation-structure formation problem. The linkage be¬ 
tween observations, especially those from the recent Cosmic Background Explorer (COBE) 
data, and theory, based on grand unified theories (GUT’s) and Friedmann-Robertson- 
Walker- (FRW-)de Sitter models, has been pursued in earnest, but most theoretical discus¬ 
sions to date are largely phenomenological and somewhat utilitarian in nature ||T4|,|T5H . This 
lack of rigor and precision is understandable for at least two reasons: the precise physical 
conditions between the Planck and GUT scales (when the most cosmologically significant 
inflationary evolutions are believed to have taken place) have not been clearly understood, 
and the theoretical framework for the treatment of processes affecting the inception and 
completion of inflation were not well developed. As stressed by one of us earlier ||Tfl| . |T7|| , the 
important physical processes which can determine whether inflation can occur, sustain, and 
finish with the necessary features are affected by at least three aspects: the geometry, topol¬ 
ogy, and dynamics of the spacetime |18|], the quantum field theory aspects pertaining to the 
analysis of infrared behavior, and the statistical mechanical aspects pertaining to nonequilib¬ 
rium processes. These quantum and statistical processes include phase transition, particle 
creation, entropy generation, fluctuation or stochastic dynamics, and structure formation 
n||2(J. Most of these invoke the quantum field and statistical mechanical aspects, and for 
processes occurring at the Planck scale (which are instrumental in starting certain models of 
inflation, such as proposed in [|6lj2l|j2^] ), also the geometry and topology of spacetime. Two 
important problems involving field theory in curved spacetime 
tion of cosmological particle creation 


18], namely, the back reac- 


!7j on the structure and dynamics of spacetimes 


2T] . [2fi| . |28| [U|. and the effects of geometry and topology of spacetime on cosmological phase 
transitions [|l7l , ^ , [36H 39ll , were investigated systematically and comprehensively in the 1970s 
and 1980s. The statistical mechanical aspect has not been considered with equal mastery. 

The statistical mechanical aspect enters into all three stages of inflationary cosmology: 
(i) At the inception: What conditions would be most conducive to starting inflation? Do 
there exist metastable states for the Higgs boson field which can generate inflation |hJD 


Can thermal or quantum fluctuations assist the inflaton in hopping or tunneling out of the 
potential barrier in the spinodal or nucleation pictures? Most depictions so far have been 
based on the finite temperature effective potential, which assumes an unrealistic equilibrium 
condition and a constant background field. However, when asking such questions in critical 
dynamics one should be using a Langevin or Fokker-Planck equation (a generalized time- 
dependent Landau-Ginzberg equation ][II[]) incorporating dynamic dissipation and intrinsic 
noise consistently, (ii) During inflation, the dynamics of the inflaton field can be more 
easily understood in terms of a Kadanoff-Migdal exponential scaling transform [52] . The 
reason why the inflaton evolves as a classical stochastic field [[43H45H at late times involves 


the process of decoherence, caused by noise and fluctuations from environmental fields K) 


this necessitates statistical mechanical considerations. The evolution of the classical density 
contrast (containing the seedings of structures) from quantum fluctuations of the inflaton 
also requires both quantum and stochastic field theory considerations |ff7| |o0|j . (iii) In the 


reheating epoch, particle creation induces dissipation of the inflaton field, and the interaction 
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of quantum fields is the source for reheating the Universe. This last epoch is the focus of 
much recent work, as we shall detail below. 

The construction of a viable theoretical framework for treating quantum statistical pro¬ 
cesses in the early Universe has been the aim of research of one of us for the past decade (for a 
review, see p7]]). This framework has now been successfully established, and its application 
to the problems mentioned above has just begun. The cornerstones are the Schwinger- 
Keldysh closed-time-path (CTP) [^-f7D| effective action and the Feynman-Vernon influence 
functional [fTlH 77[l formalisms. They are useful for treating particle creation back reaction 
|S%|,7D|1, fluctuation or noise, and dissipation or entropy problems Other essential 

ingredients include the Wigner function ||78| , |79| , the n-particle-irreducible (?rPI) effective ac¬ 
tion f|S(lEkl,l( r ~ 


, and the correlation hierarchy for treating kinetic theory processes 

and phase transition problems ||84| , |55|j . In this and three following papers we apply 
these techniques to the problems of inflaton dissipation due to parametric particle creation 
PI and reheating due to particle interaction (in the third epoch depicted above). In 


parallel, these newly developed methods in statistical held theory are now being applied to 
derive the classical stochastic dynamics of the inflaton (in the second epoch) j46|, and the 
statistical held theory of spinodal decomposition (in the hrst epoch) |4l[]. 


A. Background and issues 


All inflationary cosmologies share the feature of a period of cosmic expansion driven by a 
nearly constant vacuum energy density p (a “vacuum-dominated” era with effective equation 
of state, p = — p): In a Friedmann-Robertson-Walker (FRW) spacetime, the scale factor 
expands exponentially in cosmic time, resulting in extreme redshifting of the energy density 
of all other forms of matter and holds. As long as the interaction time scale of any physical 
process involving given helds is longer than the cosmic expansion time H -1 , the helds will 
remain in disequilibrium. This condition can prevail in all three stages of inflation, and one 
should use a fully nonequilibrium, nonperturbative treatment of the dynamics of the inflation 
held. The physics of the reheating epoch is important because it directly determines several 
important cosmological parameters which are relevant to later evolution of the Universe, 
and in principle verifiable by observational data. For example, the reheating temperature is 
a vital link between the inhationary Universe scenario and GUT scale baryogenesis 


It is generally believed that at the end of inflation, the state of the inflaton held can be 
approximately described by a condensate of zero-momentum particles undergoing coherent 
quasioscillations about the true minimum of the effective potential |[IT|JI^ , |89[| . The reheating 
problem involves describing the processes by which the many light helds coupled to the 
inflaton become populated with quanta, and eventually thermalize. It is commonly believed 
that if the helds interact sufficiently rapidly and strongly, the Universe thermalizes and turns 
into the radiation-dominated condition described by the standard Friedmann solution, but 
this has not been proven satisfactorily. 

There has been a great deal of work over the past 15 years on the reheating problem, 
and in attempting to understand reheating, a wealth of interesting physics has been revealed 
(see, e.g., plf ). To date, the work on particle production during reheating largely follows 


two distinct approaches, each pursued in two stages. 
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In the first stage of work on the reheating problem (group 1A, |9T|-|94|1 ), time-dependent 
perturbation theory was used to compute the rate of particle production into light fields 
(usually fermions) coupled to the inflaton. Particle production rates were computed in flat 
space assuming an eternally sinusoidally oscillating inflaton held. The inflaton evolution 
in FRW spacetime was modeled with a phenomenological c-number equation involving the 
Hubble parameter H and the classical inflaton amplitude cj), 


cPcj) 
dt 2 


+ m 2 (j) + (r + 3H)^- = 0, 
dt 


( 1 . 1 ) 


where T, given by the imaginary part of the self-energy of cj), is the total perturbative decay 
rate. Bose enhancement of particle production into the spatial Fourier modes of the inflaton 
fluctuation held ip (and light Bose helds coupled to the inflaton) was not taken into account. 

In the second stage of this hrst approach to the reheating problem (group IB, fl95H98[|) Eq. 
(El) was still utilized to model the mean-held dynamics, but with T computed beyond first- 


order in perturbation theory. In the work of Shtanov, Traschen, and Brandenberger |9(| and 
Kofman, Linde, and Starobinsky (KLS) |9H|, T was computed for a real self-interacting scalar 


inflaton held cj) which was both Yukawa-coupled to a spinor held -0, and bi-quadratically 
coupled to a scalar held x (KLS studied both the 0 —> —cj) symmetry-breaking and unbroken 
symmetry cases). From the one-loop equations for the quantum modes of the y, ip, and 0 
helds (in which the mean held cj) appears quadratically as an effective mass), approximate 
expressions for the growth rate of occupation numbers were derived, assuming a quasi- 
oscillatory mean held cj). For bosonic decay-product holds, it was found that hrst-order 
time-dependent perturbation theory drastically underestimates the particle production rate 
for modes which are in an instability band (due to parametric resonance). Parametric 
amplification of quantum huctuations in Bose decay-product helds can result in rapid out-of¬ 
equilibrium transfer of energy from the inflaton mean held to the (spatially) inhomogeneous 
inflaton modes and light Bose helds coupled to the inflaton. This phenomenon was called 
preheating by KLS. It has been suggested that exponential growth of quantum huctuations 
can in some cases lead to out-of-equilibrium (nonthermal) symmetry restoration in the “new” 



In both stages of this hrst approach, the back reaction of the variance of the inflaton 
on the mean-held dynamics, and of the variance on the quantum mode functions, were not 
treated self-consistently. The effect of spacetime dynamics was either excluded entirely, or 
not included self-consistently using the semiclassical Einstein equation. Due to the poten¬ 
tially large initial inflaton amplitude at the onset of reheating, particularly in the case of 
chaotic inflation |L2|, the effect of cosmic expansion on quantum particle production needs 


to be included. Since the mean held and variance (mean-squared huctuations) are coupled, 
the back reaction of particle production on the mean-held dynamics must be accounted for 
in a self-consistent manner. 

In the decade before the advent of inhationary cosmology, there was active research on 
quantum processes in curved spacetimes. An important class of problems is vacuum particle 
creation K3Lp7| and its effect on the dynamics and structure of the early Llniverse 
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at the Planck time. The effect of spacetime dynamics and the importance of parametric 
amplification on cosmological particle creation were realized very early |23| 


27 . Most of 


the effort in the latter part of the 1970s was focused on obtaining both a regularized energy- 
momentum tensor and a viable formalism for the treatment of back reaction effects. The 
wisdom gained from work in that period before the inflationary cosmology program was 
initiated is particularly relevant to the reheating problem. Simply put, for obtaining a finite 
energy-momentum tensor for a quantum field in a cosmological spacetime, the adiabatic 


27,102-104] and dimensional 


formalism |58|j60 
method 


105|| regularization methods are the most useful. For studying 

effective action 
in-out”) 


the back reaction of particle creation, the Schwinger-Keldysh (CTP, “in-in 

a is more appropriate than the usual Schwinger-DeWitt 


IoWotT . 


The second approach to the post-inflationary reheating problem is built upon the body 
of earlier work on cosmological particle creation. Following the application of closed-time- 

several authors 


with cubic 


coupled to fermions 1109 


path techniques to nonequilibrium relativistic field theory problems 
(which we call group 2A) derived perturbative mean-field equations for a scalar inflaton 
TAJ and quartic ||108|| self-couplings, as well as for a scalar inflaton Yukawa- 
The closed-time-path method yields a real and causal mean-held 
equation with back reaction from quantum particle creation taken into account. For the 
case of Bose particle production, perturbation theory in the coupling constant is known to 
break down for sufficiently large occupation numbers, which occurs on the time scale T\ for 

mm 


parametric resonance effects to become important 


It is, therefore, necessary to 


employ nonperturbative techniques in order to study reheating in most inhationary models. 

The second stage of work in this second approach to the reheating problem used the 
closed-time-path method to derive self-consistent mean-held equations for an inflaton cou¬ 
pled to lighter quantum helds (group 2B, ||112| - |118| ]). In the hrst of these studies ||1 1 2| |115|| , 
the coupled one-loop mean-held and mode-function equations were solved numerically in 
Minkowski space, implicitly carrying out an ad hoc nonperturbative resummation in h. In 
the one-loop equations, the variances for the inflaton (<p 2 ) and light Bose helds (y 2 ) do 
not back-react on the mode functions directly. However, mean-held equations were derived 
for an 0(iV)-invariant linear a model (with a A$ 4 self-interaction) at leading order in the 
large -N approximation by Boyanovsky et al. ||110|| . In this approximation, the variance does 


back-react on the quantum mode functions. At leading order in the 1/N expansion, the 
unbroken symmetry dynamical equations for the quartic O (N) model are formally similar 
to the dynamical equations for a single A<h 4 held theory in the time-dependent Hartree-Fock 
approximation [^Cj. The nonequilibrium dynamics of the quartically self-interacting O (N) 


held theory in Minkowski space has been numerically studied at leading order in the 1/N ex- 

mmm 


pansion in both the unbroken symmetry 
cases. Some analytic work has been done on the self-consistent Hartree-Fock mean-held 


equations for a quartic scalar held in Minkowski space [111]. In addition, the Hartree-Fock 


equations for a Ad* 4 held in the slow-roll regime have been studied numerically in Minkowski 

ra 


space 


and in FRW spacetime 


However, the effect of spacetime dynamics on re¬ 


heating in the 0(A r ) held theory has not (to our knowledge) been studied using the coupled, 
self-consistent semiclassical Einstein equation and matter-held dynamical equations, though 
some simple analytic work has been done on curvature effects in reheating |p7| , |116 


semiclassical equations for one-loop reheating in FRW spacetime were derived in ||117 


The 

The 
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( f> 2 x 2 theory has been studied in FRW spacetime by [|54| , |1181 , |123|| . In addition, numerical 
work has been done on symmetry-breaking phase transitions in both a A$ 4 scalar held in 
de Sitter spacetime [ |124|| , and an 0(7V) theory in FRW spacetime ||125| , |126|| . The group 2B 
studies represent the state of the art on the preheating problem. 


B. Our work how it differs from others 


In the present work we study the nonperturbative, out-of-equilibrium dynamics of a 
minimally coupled scalar 0(N ) held theory (with quartic self-interaction) in a spatially hat 
FRW universe whose dynamics is given self-consistently by the semiclassical Einstein equa¬ 
tion. The purpose of this study is to understand the preheating period in inhationary cos¬ 
mology, with particular emphasis on the effect of spacetime dynamics on the phenomenon 
of particle production via parametric amplification of quantum huctuations. Of primary 
interest is obtaining the dynamics of the inhaton (including back reaction from created 
particles) using rigorous methods of nonequilibrium held theory in curved spacetime [|69 


We have chosen to focus in this paper on parametric amplification of quantum huctuations 
because this phenomenon can be the dominant effect in the preheating stage of unbroken 
symmetry inhationary scenarios, among which the chaotic inflation scenarios most directly 
necessitate (through initial conditions) considerations of Planck-scale physics. “New” inha¬ 
tionary scenarios which involve a spontaneously broken symmetry often contain additional 
subtleties (e.g., infrared divergences, spinodal instabilities), and are the subject of ongoing 
investigation : |4 T} F(j| . The results of our work are, therefore, particularly relevant to chaotic 


inflation scenarios |127 [. The additional interactions which should be included to treat the 
broken-symmetry case are discussed in [l|. 

In the preceding paper []I|J we have derived the evolution equations for the mean held 
($ H ) (subscript H denotes the Heisenberg held operator) and mean-squared huctuations 
(variance) ($^) — ($ H ) 2 using the closed-time-path (CTP), two-particle-irreducible (2PI) 
effective action [Q in a fully covariant form. Here we use these results for the case of 
spatially hat FRW spacetime. The quantum state for the held theory (in the case of FRW 
spacetime) consists of a coherent state for the spatially homogeneous held mode, and the 
adiabatic vacuum state for the spatially inhomogeneous modes. At conformal past infinity, 
the spacetime is assumed to be asymptotically de Sitter, and the mean held is chosen to be 
asymptotically constant. 

In this paper we study the O(N) held theory using the 1/N expansion, which yields 
nonperturbative dynamics in the regime of strong mean held. This is particularly important 
for chaotic inflation scenarios [jfj , in which the inhaton mean-held amplitude can be as large 
as Mp/3 at the end of the slow-roll period ||T2| , |128|| . Treatments of the reheating problem 
which rely on time-dependent perturbation theory do not apply to such cases where the 
inhaton undergoes large-amplitude oscillations, in contrast with nonperturbative methods 
such as large N. 

We now summarize the principal distinctions between our work and previous treatments 
of preheating in inhationary cosmology. Our work improves on the group 1A methodology 
by including parametric resonance effects. As it is based on hrst-order, time-dependent 
perturbation theory, the group 1A approach cannot correctly describe the inhaton dynamics 
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with large initial mean-field amplitude. In addition, our work improves on both the group 1A 
and group IB studies by including the effect of back reaction from quantum particle creation 
on both the mean field and the inhomogeneous modes. We are treating the inflaton dynamics 
from first principles, without assuming a phenomenological equation (with a damping term 
T(f> put in by hand) for the mean field. In our approach the damping of the mean field is due 
to back reaction from quantum particle production in the self-consistent equations for the 
mean field and its variance. In contrast, the analytic results of the group 1A and IB work 
are based on the assumption of either large-amplitude mean-field oscillations (A0 2 /2 m 2 ) 
or harmonic oscillations (m 2 3> </> 2 /2), and, therefore, cannot describe the interesting case 
of inflaton dynamics in which neither term dominates the tree-level effective mass, i.e., 
m 2 ~ \(j) 2 . Furthermore, our work improves on group 1A, group IB, and group 2A in that 
the closed-time-path effective action is computed in curved spacetime without assuming that 
H tq (where To is the period of mean-field oscillations). In our work, the dynamics of 
the two-point function (which reflects quantum particle production) is formulated in curved 
spacetime assuming only that semiclassical gravity is valid, i.e., M P tU H. 

Most significantly, our work improves on all the previous treatments in that it includes 
curved spacetime effects systematically using the coupled, self-consistent semiclassical Ein¬ 
stein equation and matter field equations. Among the group 2B studies of preheating 
dynamics, inflaton dynamics has been studied primarily in fixed background spacetimes: 
Minkowski space ||1Q1|J1 10| 


de Sitter space [|124|J126|Jl29|l , and in radiation-dominated, spa¬ 


tially flat FRW spacetime 1126]. In the present work, the spacetime is dynamical, with 


the renormalized trace of the semiclassical Einstein equation governing the dynamics of the 
scale factor a. This permits quantitative study of the transition of the spacetime from the 
(slow-roll) period of vacuum-dominated expansion to the radiation-dominated (“standard”) 
FRW cosmology. In particular, our method yields the spacetime dynamics naturally, with¬ 
out making reference to an “effective Hubble constant” (which has been used in calculations 
on a fixed background spacetime [|129|| ). 

With additional couplings (see [(TJ), our method may also be used to study preheating in 
“new” inflationary scenarios ||. In new inflation, the vacuum-dominated expansion of the 
Universe is typically driven by the classical potential energy of the mean field as it rolls to¬ 
wards the symmetry-broken ground state. In one of the group 2B studies (Boyanovsky et al. 

129||), a quench-induced phase transition is studied with small initial mean-field amplitude, 


in which the classical terms in the mean-field equation are dominated by spinodal fluctua¬ 
tions. As a result, the mean field in their model does not oscillate about the symmetry-broken 
ground state as is generally expected in a new inflation preheating scenario (this point was 


emphasized in 


The initial conditions studied in [129|| are more appropriate to a study 


of defect formation in a quench-induced phase transition than preheating dynamics in new 
inflation. 

is not generally covariant (as 


In addition, the renormalization scheme employed in [129 


can be seen by comparing it with ||131|| ), and covariant conservation of the renormalized 
energy-momentum tensor is put in by hand. The regularization scheme employed here is 
the well-tested adiabatic regularization lp7| .|102| [109]] , which is simple to use and physically 
intuitive. It also ensures both covariant conservation of the regularized energy-momentum 
tensor and agreement with manifestly covariant regularization procedures such as point 
splitting ||131| . 















































A related difference between our approach and that of some of the group 2B studies is 
the choice of vacuum state. The choice of initial conditions for the quantum mode functions 
in most studies of reheating in FRW spacetime ||118| , |122| , |132|| has been to instantaneously 
diagonalize the matter-held Hamiltonian at the initial-data hypersurface. However, as has 
been pointed out long ago ||133|| , this method does not correspond to the vacuum state which 
registers the least particle flux on a comoving detector. In our work we use the de Sitter- 
invariant (or Bunch-Davies) vacuum, obtained via the adiabatic construction; the adiabatic 
vacuum most closely aligns with an intuitive notion of vacuum state in a cosmological 
spacetime |TE| . 


In many of the group 2B studies ||110|J126| , |129|| the large -N equations for the mean held 
and variance are derived using a factorization method which does not readily generalize to 
next-to-leading order in the 1/N expansion. In all of the group 2B studies of nonperturbative 
inhaton dynamics of which we are aware, the equations for the mean held and variance are 
not derived using methods which encompass higher-order correlations in the Schwinger- 

mm 


Dyson hierarchy. As found in earlier studies of phase transitions 


this is necessary 


in order to derive the correct infrared behavior of a quantum held in a study of critical 
phenomena. In the present work, we use the result of the preceding paper |3J in which the 
CTP two-particle-irreducible (2PI) formalism is derived. It has a direct generalization in 
terms of the n-particle-irreducible (nPI) “master” effective action j82|. The master effective 


action can be used to derive a self-consistent truncation of the Schwinger-Dyson equations 

HI 


to arbitrary order in the correlation hierarchy 


The techniques employed here are, 


therefore, most readily generalized to the study of phase transitions in curved spacetime, 
where higher-order correlation functions can become important [Pi] . 

In summary, our approach to the inhaton dynamics problem has the following advantages: 
it is nonperturbative and fully covariant; it is based on rigorous methods of nonequilibrium 
held theory in curved spacetime; we use the correct adiabatic vacuum construction; and we 
employ an approximation scheme which can be systematically generalized beyond leading 
order, within a fully covariant and self-consistent theoretical framework. 


C. Summary of major findings 

Our results are obtained by solving the mean-held and spacetime dynamics self- 
consistently using the coupled matter-held and semiclassical Einstein equations in a FRW 
spacetime, including the effect of back reaction of the variance on the mean held. Within the 
leading-order, large N approximation used here, we find that (using the conventional value 
for the self-coupling, A = ICE 14 ) for sufficiently large initial mean-held amplitude, paramet¬ 
ric amplification of quantum fluctuations is not an efficient mechanism of energy transfer 
from the mean held to the inhomogeneous held modes. In this case the energy density of 
the inhomogeneous modes remains negligible in comparison to the mean-held energy density 
for all times. This can be understood from the time scales for the competing processes of 
parametric resonance and cosmic expansion. When the time scale for parametric amplifica¬ 
tion of quantum fluctuations T\ is of the same order as (or greater than) the time scale for 
cosmic expansion H -1 , cosmic expansion redshifts the energy density of the inhomogeneous 
modes faster than it increases due to parametric resonance. We find that this occurs when 
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(j) > Mp/300, for the model and coupling studied here. 

In many chaotic inflation scenarios, the mean-held amplitude at the end of the slow- 
roll period can be as large as M P /3 Hl2l,|128|| . In light of our result, in such models, it 


is clearly essential to include the effect of spacetime dynamics in order to study mean-held 
dynamics and resonant particle production during reheating. In addition, our result indicates 
that for the case of a minimally coupled A<h 4 inhaton with unbroken symmetry, parametric 
amplification of its own quantum fluctuations is not a viable mechanism for reheating, unless 
the coupling is significantly strengthened (see |55| , |56| ). Parametric amplification of quantum 
fluctuations may still play a dominant role in the reheating of chaotic inhaton models with 
an inhaton coupled to other helds, e.g., a 4> 2 x 2 model. This is the subject of a forthcoming 
paper |86 . 

For more moderate cosmic expansion, where H~ l > 100 Ti, parametric amplification 
of quantum fluctuations is an efficient mechanism of energy transfer to the inhomogeneous 
modes, and the asymptotic effective equation of state is found to agree with the prediction 
of a two-huid model consisting of the elliptically oscillating mean held and relativistic energy 
density contained in the inhomogeneous mode occupations. In a collisionless approximation, 
the mean held eventually decouples from the mean-squared huctuations (variance) and at 
late times undergoes asymptotic oscillations which are damped solely by cosmic expansion 
|128|| . For the case when cosmic expansion is subdominant, H -1 7i, the mean-held 

dynamics and the growth of quantum huctuations are in agreement with results of studies 
of preheating in Minkowski space [|101|| . In particular, the total adiabatically regularized 
energy density is found to be constant (to within the limits of numerical precision) for the 
case of H > oo, in agreement with the predictions of held theory in Minkowski space. 

While there has been a large volume of work on the preheating period and particle pro¬ 
duction, the thermalization of inhationary models has not yet been understood from hrst 
principles. From the nonperturbative inhaton dynamics, Boyanovsky et al. claimed that 
H3 the Boltzmann equation is inadequate for studying collisional thermalization at the 
end of the preheating stage. In particular, in the leading-order, large -N approximation em¬ 
ployed here (and in the one-loop approximation which it contains), this model also does not 
thermalize. However, it still may approach a radiation-dominated effective equation of state 
in Minkowski space as found in ||101|| . Clearly, a hrst-principles analysis of thermalization 
is necessary. Continuing the early work of kinetic held theory J69|], and the recent work on 
correlation hierarchy [^], we know that such a hrst-principles analysis should involve at a 
minimum the full two-loop, two-particle-irreducible effective action (or alternatively, next- 
to-leading order in the large -N approximation). Since it represents a rigorous truncation 
of the full Schwinger-Dyson hierarchy, in this sense it is the most natural generalization 
of the collisionless approximations used previously to study reheating. However, the equa¬ 
tions derived from it for the mean-held and gap equation are nonlocal and hence difficult to 
solve even numerically ||119|1 . This paper is, therefore, concerned only with preheating via 


parametric resonance particle creation; work on thermalization is in progress |87| . 
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D. Organization and notation 


This paper is organized as follows. In Sec. |TI] we present the general theory of noneqni- 
librium dynamics of a scalar field in curved spacetime, including a summary discussion of 
reheating in inflationary cosmology. In Sec. 0 we specialize to the case of spatially flat 
FRW spacetime, and derive the dynamical equations. The results of numerically solving the 
evolution equations are contained in Sec. |IV C| . Discussion and conclusions follow in Sec. 0. 

Throughout this paper we use units in which c = 1. Planck’s constant h is shown 
explicitly (i.e., not set equal to 1) except in those sections where noted. In these units, 


Newton’s constant is G = HM P 2 , where M P is the Planck mass. We work with a four¬ 
dimensional spacetime manifold, and follow the sign conventions]]] of Birrcll and Davies (18 


for the metric tensor g fiU , the Riemann curvature tensor R lwap , and the Einstein tensor G pv . 
We use greek letters to denote spacetime indices. The beginning latin letters a, b, c, d, e, / 
indicate the time branch (see [I|, Sec. II), and the middle latin letters i,j, n are 

reserved as indices in the 0(1V) space (see Sec. ED - Einstein summation convention over 
repeated indices is employed. 


II. A<h 4 INFLATON DYNAMICS IN FRW SPACETIME 

A. A<f > 4 quantum fields in curved spacetime 

As a simple model of inflation, let us consider a scalar A® 4 field in semiclassical gravity, 
where the matter field is quantized on a classical, dynamical background spacetime. The 
classical action has the form 


s[<f>,gn = F\ s r\+s*[<i>, s r\i 

where S F is the matter field action, which for the scalar A® 4 theory takes the form 


S F [<P,g > “'] = -- / d'xsf^j 


<p(n + m 2 + £R)<I> + -W 


( 2 , 1 ) 


( 2 , 2 ) 


and for a renormalizable theory, the gravity action S G must have the form (18( 1135 
1 


s c \ g n = 


167tG 


d 'xy/^j [R- 2A + cR 2 + bR a(3 R a p + aR a ^ s R aPlS \ . (2.3) 


In Eqs. ( p . 2| ) and Q2.3|) , A is the coupling constant (with dimensions of inverse mass times 
inverse length), m is the “mass” (with dimensions of inverse length), £ is the dimensionless 
coupling to gravity, G is Newton’s constant (with dimensions of length divided by mass), 
and a, b, and c are constants with dimensions of length squared. The symbol □ denotes the 


1 In the classification scheme of Misner, Thorne, and Wheeler [134], the sign convention of Birrell 
and Davies [18] is classified as (+,+,-!-). 


11 
























Laplace-Beltrami operator in terms of the covariant derivative V M , R is the scalar curvature, 
R^ v is the Ricci tensor, R a 0 ^s is the Riemann tensor, and A J—g is the square root of minus 
the determinant of the metric tensor g^ v . 

The inflaton held 0 is then quantized on the classical background spacetime; we denote 
the Heisenberg held operator by <f> H , and the quantum state by 10). Of particular interest 
in a study of inhaton dynamics are the mean held 


00) = (0|$hO)|0>, 

(2.4) 

the huctuation held 


<PhO) = $h0) - 00)) 

(2.5) 

and the mean-squared fluctuations, or variance 


(0OhO)I0) = (0I$hO)I0) - OOhOM 2 . 

(2.6) 


In a previous paper |T|, a systematic procedure was presented for deriving real and causal 
evolution equations for the mean held, two-point function, and the metric tensor in semiclas- 
sical gravity. Assuming a globally hyperbolic spacetime, one can evolve the coupled evolution 
equations forward from initial data specihed at an initial Cauchy hypersurface. The evolu¬ 
tion equations follow from functional differentiation (and subsequent held identifications) of 
the closed-time-path (CTP) two-particle-irreduciblc (2PI) effective action, r[0 ± , G±±, g±]. 
The CTP-2PI effective action is a functional of the mean held 0, two-point function G , and 
metric tensor g lllJ , which now carry not only spacetime labels but also time branch labels, 
which have an index set {+,—}• The evolution equations for 0, {(p 2 ), and g^ then follow 
from 


5(5 G bf]+r[0 ± ,G ±± ,^f]) 


5r [0 ± ,G ±± ,<7f] 


Ha 

ST[4> ± 1 G±±,g^} 


SG, 


ab 



= 0, 

(2.7a) 

i + =i_=P, gf =g* u =g^ 




= 0, 

(2.7b) 

<A+= 0 -=<A; g 1 ^ =g^ =g> l ' J 




= 0. 

(2.7c) 

i+=i-=t; gf =g^ =g^ 




The variance is related to the coincidence limit of CTP two-point function as follows: 


(P(x)) = hG a b(x,x) 


( 2 , 8 ) 


for all a, b G {+, —The energy-momentum tensor (T0,) is dehned by 


(r„„) = 


<sr[0±,G ±± , 9 £ 


(IIP I 




5g 


fLV 


(2,9) 


$ + =4,_=$-, g^ , =gt v =g» v 


which (after renormalization) enters as the source of the semiclassical Einstein held equation: 
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87rG(T^). 


( 2 . 10 ) 


GW = 


Eqs. (|2.7a|) - (|2.7cj ) constitute a set of coupled, nonlocal, nonlinear equations for the mean 
field, two-point function, and metric tensor. The renormalized versions are what enter into 
the description of inflaton dynamics. The CTP-2PI effective action can be computed using 
diagrammatic methods described in the previous paper, where a covariant expression for T 
was computed in a general curved spacetime (truncated at two loops). 


B. Inflaton dynamics in FRW spacetime 

We now consider a spatially flat Friedmann-Robertson-Walker (FRW) spacetime, which 
is spatially homogeneous, isotropic, and conformally flat. Its line element can be written in 
the form 


ds 2 


a (v) 2 


3 

drf — ^(cfc*) 2 


2=1 


( 2 . 11 ) 


where a is the scale factor, x l (i e {1,2,3}) are the physical position coordinates on the 
spatial hypersurfaces of constant conformal time rj (related to the cosmic time t by r] = 
f dt/a ). The Hubble parameter, which measures the rate of cosmic expansion, is 


H{v) = ( 2 . 12 ) 

a 

where the over-dot denotes differentiation with respect to cosmic time t. Given our choice 
of sign convention and metric signature, the Ricci tensor in the FRW coordinates is given 
by 


Roo — 3 


Rij — 


a 

a 


,A2 


a ( a ' 

— + 


A2 


Oij, 


(2.13a) 

(2.13b) 


where the prime denotes differentiation with respect to r), and i?oo is the component of the 
Ricci tensor proportional to dr] <g> dr]. The scalar curvature is 


and the Einstein tensor is 


R 



Goo — 
Gij — 


3(a') 2 
a 2 ’ 

'2a" {a') 2 ' 

a a 2 


jij- 


Finally, the volume form on M is 


(2.14) 


(2.15a) 

(2.15b) 
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e M = a 4 (dr] A dx 1 A dx 2 A dx 3 ). (2-16) 

The higher-order (e.g., R 2 ) geometric terms in the geometrodynamical field equation are not 
shown because the renormalized constants b and c are set to zero in Sec. pTTDj. 

In restricting the spacetime to be a spatially flat FRW, we are reducing the number of 
degrees of freedom in the metric: 


9 tLv -> o(r 


(2.17) 


This reduction should not be carried out in the 2PI generating functional T[0, G, g^], but 
only in the equations of motion (|2 . 7a|) — (|2 . 7c|) . This is because functional differentiation of 
T [</>, G, a~ 2 ri^ v ] with respect to the scale factor a gives only the trace of the energy-momentum 
tensor, a~ 2r q^ u {T^ v ), and not the additional constraint equation which the initial data must 
satisfy. 

The spatial homogeneity and isotropy of FRW spacetime permits only two algebraically 
independent components of the energy-moment urn tensor, which in the FRW coordinates 
of Eq. ([2. 11| ) are given by (T 00 ) and (T.n)] all other components are zero. These must be 
functions of rj only (due to spatial homogeneity). For the purpose of numerically solving the 
semiclassical Einstein equation, it is convenient to work with the trace 


T = 


(2.18) 


instead of (T tl ). The trace T enters into the dynamical equation for a(rj), and (T 00 ) enters 
into the constraint equation. 

Another consequence of the spatial symmetries of FRW spacetime is the restriction on 
the generality with which we may specify initial data for dynamical evolution. Let us choose 
to specify initial data on a Cauchy hypersurface E r)0 of constant conformal time t] 0 . In the 
Heisenberg picture, 0 for consistency with spatial homogeneity, the quantum state \(j>) must 
satisfy 

(0|$ H (^o,^)|0) = 0M, (2.19a) 

((p\K(Vo,x)\(p) = 0'(?/o), (2.19b) 

for all x G M 3 , where <F H is the Heisenberg field operator for the scalar field. The values of 
<p(r]o) and (A (bo) constitute initial data for the mean field. In addition, the quantum state 
must satisfy 


(0|( / 9 h (??o,^)Ah(^o,^ , )|0) = F(Vo, l^-^l), (2.20a) 

vp = F'(r) 0 , \x - f'|), (2.20b) 

dr 1 \vo 


in terms of an equal-time correlation function F(rj o, \x — T'|) which is invariant under si¬ 
multaneous translations and rotations of x and x'. As defined in Eq. (j2.5p, <p H denotes the 


2 As discussed in Sec. IIC , for our purposes it is sufficient to consider only the case of a pure state. 
The analysis can, however, be easily extended to encompass a mixed state with density matrix p. 
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Heisenberg field operator for the fluctuation field. The spatial Fourier transform of F is 
related to the power spectrum of quantum fluctuations at rj 0 for the quantum state |0). 
Alternatively, we may say that F(r] 0 ,r) and F'(r] 0 ,r ) give initial data for the evolution of 
the two-point function G ++ via the gap equation (|2.7c|) . The symmetry conditions fl2.19a|) , 
( j2.19b|) , ( f2.20a| ), (|2.20b|) , along with the spatial symmetries of the classical action in FRW 
spacetime, guarantee that the mean field and two-point function satisfy spatial homogeneity 
and isotropy for all time, i.e., 


( < h H (^)) = to), 

G++(x, x') = G ++ (pk ?/, \x - x'\), 


(2.21a) 

(2.21b) 


for all x G M. The conditions (|2.21a|) , (|2.21b|) permit a formal solution of the gap equation 
( p.7c|) for G ++ in terms of homogeneous mode functions, via a Fourier transform in comoving 


momentum k , as shown in Sec. |I1I B| . By rotational invariance, the Fourier transform depends 

only on the magnitude k = \/k ■ k. Of course, the quantum state \4>) is not uniquely defined 
by the spatial symmetries; a unique choice of the initial conditions for <p and G a b at E 


-“no 


is 


(in the Gaussian wave-functional approximation) equivalent to choosing 
quantum state depends on the physics of the problem we wish to study. 

As a consequence of covariant conservation of the energy-momentum tensor 

V^(T^) = 0, 

the functions (T 0 o(?7)) and ('2~h(r/)) satisfy 

dr ,\ (Ta) d 


The choice of 




a 2 dr] 


( 2 . 22 ) 


(2.23) 


which comes from taking the v = 0 component of Eq. ( p.22| ) . I 11 analogy with the continuity 
relation for a classical perfect fluid in FRW spacetime, 


Gd - ~ 4 , G) • 

we may define the energy density p and pressure p, by 


(2.24) 


P(v) = — ( T oo(v)), 

cr 

p 4) = -ATM)- 


(2.25a) 

(2.25b) 


However, the quantity p should not be interpreted as the true hydrodynamic pressure until 
a perfect-fluid condition is shown to exist; otherwise, bulk viscosity corrections can enter 
into Eq. ( |2.25b|) ||136|| . The effective equation of state is defined as a time average (over the 


time scale T\ for the matter field dynamics, to be discussed in Sec. IID ) of the ratio p/p, 

(2.26) 


_ = P 
~ P 
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The effective equation of state 7 (where the bar denotes a time average) is an important 
quantity in differentiating between the various stages of inflationary cosmology. 

Several solutions to the semiclassical Einstein equation ( [2.10| ) for idealized equations 
of state are of particular interest in cosmology. The effective equation of state 7 = — 1 
(eternally “vacuum dominated”) leads to a solution a(i]) = —1 /(Hi]), for — 00 < r) < 0, 
where H = 8nGp/3 and p is a constant. This solution corresponds to the “steady-state” 
coorclinatization covering one-half of the de Sitter manifold |TS|. The effective equation of 
state 7 = 0 corresponds to nonrelativistic matter, in which case the scale factor conformal¬ 
time dependence is a oc if. The effective equation of state 7 = 1/3 corresponds to relativistic 
matter, and its scale factor conformal-time dependence is a oc 7 . 


C. Initial conditions for post-inflaton dynamics 


In most realizations of inflationary cosmology, the Universe evolves through a period 
in which a dominant portion of the energy density p comes from a quantum field $ H , the 
inflaton field, whose effective equation of state [defined as in Eq. ( |2.26| )1 is 7 ~ — 1. In 
chaotic inflation, this condition is due to the fact that the inflaton held is in a quantum 
state |0) in which the Heisenberg held operator <f> H acquires a large (approximately spatially 
homogeneous) expectation value, dehned by 


4>(v) = (0|$ H (a;)|0)- 


(2.27) 


A requirement for chaotic inflation is that the potential energy U(0) of the expectation 
value 0 dominates over both the spatial gradient energy [coming from ((Vy? H ) 2 )] and kinetic 
energy for the inflaton held, and the energy density of all other quantum helds coupled 
to the inflaton. The potential energy U(0) gives a contribution to the energy-momentum 
tensor satisfying precisely 7 = —1. During inflation, the scale factor grows by a factor 
of approximately exp (//At), where At is the interval of inflation in cosmic time, typically 
larger than 60// . While the Universe is inflating, the expectation value (<I> H ) is slowly 

rolling toward the true minimum of the effective potential. (In reality, the situation is much 
more complicated than this. The effective potential is an inadequate tool for studying out- 
of-equilibrium mean-held dynamics |T6| , |137|| .) Assuming the Universe was in local thermal 
equilibrium prior to inflation, the temperature during inflation decreases in proportion to 
1/a. The energy density of any relativistic (nonrelativistic) helds coupled to the inflaton is 
proportional to 1/a 4 (1/a 3 ). The contribution to the quantum energy density from spatial 
gradients of fluctuations about the inflaton held is proportional to 1/a 4 (see Sec. P~TT| below). 
Most importantly, any inhomogeneous modes Sfik of the mean field which might exist at the 
onset of inflation are redshifted. The physical momentum of a quantum mode, k phyB = k/a, 
decreases as 1/a relative to the comoving momentum k. The quantum state of any held 
coupled to the inflaton at the end of inflation is, therefore, approximately given by the 
vacuum state. The inflaton held is well approximated by a spatially homogeneous mean 
held, with vacuum fluctuations around the mean-held configuration. The mean held can 
be thought of as representing the coherent oscillations of a condensate of zero-momentum 
inflaton particles. 
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Let us consider the case of inflation driven by a single self-interacting scalar field 0 (with 
unbroken symmetry) in spatially flat FRW spacetime. The above arguments imply that one 
can model post-inflationary physics with a quantum state |0) which at 7] 0 corresponds to a 
coherent state for the held operator <f> H [in which ( 0 |<f> H (x)| 0 ) = and the fluctuation 

held (/? H = <f> H —0 is very nearly in the vacuum state.f] Then for rj < r ] 0 , (T 00 ) is dominated by 
the classical energy density of the mean held 0. The 00 component of the Einstein equation 
then yields 


a' 18nGp c 

a* = V 3 ’ 

where p c is the classical energy density of the mean held, defined by 

Pc = (i'f + vm 


The mean held 0 satishes the classical equation 

0" + —0' + aV'(0) = 0, 

a 


(2.28) 


(2.29) 


(2.30) 


where U(0) denotes the classical potential. For the A<f > 4 theory, the potential is [from the 
Minkowski-space limit of Eq. ( [2.11) 1 

E(0) = I m 2 0 2 +^0 4 . (2.31) 

The assumption that the Universe is inflating (i.e., 7 ~ —1) for rj < r/o requires that the 
energy density p c be potential dominated, 

V{i) » (2.32) 

and that the mean held satishes the slow-roll condition, 



a 


(2.33) 


Given Eqs. (|2.32|) and 


an approximate “0th adiabatic order” solution to the Einstein 


equation can be obtained [normalized to 0 ( 70 ) = 1 ], 


a ^ 1 + 


(2.34) 


3 Though this is a pure quantum state, the methods employed in this study can be used to treat a 
quantum field theory in a mixed state (for example, a system initially in thermal equilibrium with 
a heat bath). 
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where H is a slowly varying function of 77, given by 


H(v) = 


8tt Gp c (i]) 


(2.35) 


From Eq. (|2.35|) , we can evaluate the expansion rate nonadiabaticity parameter [@ for 

). During slow-roll it follows from conditions (|2.32|) and ([2.33D that 


77 < 770 using Eq. 


».-s= 




< 1 . 


32?rG 


v(j>y 


(2.36) 


The solution ( |2.34|) for 0(77) is exact in the limit of constant H (corresponding to a constant 
0 at the tree level). For simplicity, let us assume that 0 goes to a constant value > M P 
in the asymptotic past, 77 —> — 00 . The spacetime is then asymptotically de Sitter, with 
the scale factor having an asymptotic cosmic-time dependence a(t) — exp (Ht). Because 
the enormous cosmic expansion during the slow-roll period redshifts away all nonvacuum 
energy in the Universe, it is reasonable to assume that the quantum state |0) would register 
no particles for a comoving detector coupled to the fluctuation held < 7 ? at conformal-past 
infinity; i.e., that the fluctuation held ip is in the vacuum state at 77 —> — 00 . This would 
mean that a ~ l/(Hr]) at 77 —> — 00 . This spacetime is not asymptotically static in the 
past, but it is conformally static with a conformal factor whose nonadiabaticity parameter 
vanishes at conformal-past infinity. Therefore, the best approximation to a “no-particle” 
state for a comoving detector in the asymptotic past is given by the adiabatic vacuum || 102 | 
constructed via matching at 77 —> — 00 . All higher-order adiabatic vacua will in this case 
agree at conformal past infinity. 

To construct the 7 ith-order adiabatic vacuum matched at an equal-time hypersurface 
one hrst exactly solves the conformal-mode function equation for the quantum held 




[see Eq. ( |3.7| ) below]. Since the mode-function equation is second order, each k mode has 
two independent solutions, which can be represented as Uk and u k . A particular solution 
consists of a linear combination of Uk and u k . The adiabatic vacuum is constructed by 
choosing (for each k ) a linear combination which smoothly matches the 77 th-order positive 
frequency WKB mode function at . The resulting orthonormal basis of mode functions 
is then used to expand the Heisenberg held operator <h H (a;) in terms of a*, and a\. The 
vacuum state is dehned by a 0 vac) = 0 for all k , which can be shown to correspond (in the 
Vm —> —00 limit) to the de Sitter-invariant, adiabatic (Bunch-Davies) vacuum. 


D. Post-inflation preheating 


Inflation ends when the mean held has rolled down to the point where condition ( j2.33| ) 
ceases to be valid, which we assume occurs at conformal time r] 0 . The inflaton mean held then 
begins to oscillate about the true minimum of the effective potential, leading to a change in 
the effective equation of state. A harmonically oscillating scalar mean held (m 2 ^ A0 2 /6) 
has an effective equation of state 7 = 0 , and a scalar inflaton undergoing extreme large- 
amplitude oscillations (A0 2 /6 m 2 ) has an effective equation of state 7 = 1/3 ]Tl|. In 
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realistic models, the inflaton field is coupled to various lighter fields consisting of fermions 
and/or bosons. These couplings, as well as the inflaton’s self-coupling, provide mechanisms 
for damping of the mean-field oscillations via back reaction from quantum particle produc¬ 
tion, and energy transfer to the lighter fields and the inflaton’s inhomogeneous modes. 

Let us consider the scalar A<f> 4 field theory with unbroken symmetry in Minkowski space 
[with classical action given by the Minkowski-space limit of Eq. ( |2.1|) 1, and suppose that 
the mean field 0 oscillates about the stable equilibrium configuration 0 = 0 with initial 
amplitude 0 O . For the moment we are neglecting the effect of spacetime dynamics, i.e., 
assuming a{jj) = 1. The time scale for oscillations of the mean field is given by m 


T 0 = 


4 K(k) 

iVT+P 


(2.37) 


where / and k are defined by 


/ = \ 

o m 


k = 


f 


VW+P)' 


(2.38a) 

(2.38b) 


and K(k ) is the complete elliptic integral of the first kind [|138|| . For harmonic oscillations 
where 


^0o < 


(2.39) 


time-dependent perturbation theory was used in the group 1A studies (see Sec. |) to com¬ 
pute the damping rate T for the mean field in the A<f> 4 model. At lowest order in A, the 
damping rate for the mean field 0 corresponds to the rate for four zero-momentum, free-field 
excitations of the inflaton to decay into a ip (fluctuation field) particle pair, due to the A0 4 
self-coupling | |)6|| , 


r* ^ o(i) 


Wo ) 2 

47rm 


(2.40) 


with vacuum initial state for <p. The symbol 0(1) denotes a constant of order unity. In 
addition to the assumption ( |2.39| ), there is another crucial assumption in the derivation of 
Eq. ([2.40| ), namely, that the dominant contribution to the decay rate is given by the lowest 
order, |vac) —> |1^, 1_^) process, where the occupation numbers are for the fluctuation field 


<p. ft can be shown [|23| , |139|| that for this (bosonic) case, the perturbative decay rate into the 
k momentum shell for the fluctuation field ip is enhanced by 1 + 2 n when the occupation of 
the k shell is n. This is a stimulated emission effect due to Bose statistics.^ The use of Eq. 


1 In contrast with the case with Bose fields, the use of time-dependent perturbation theory to 
study inflaton decay into fermions via a Yukawa coupling does not require the condition nr -C 1, 
because of the Pauli exclusion principle. It is still necessary, however, to assume weak coupling (or 
small mean-field amplitude) in order to use perturbation theory |{94||)G[|. 
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( |2.4C| ) to estimate the damping rate thus implicitly assumes that for all k, the fluctuation held 
occupation numbers are small, i.e., n% -C 1. This is because time-dependent perturbation 
theory in terms of the A$ 4 interaction corresponds to an expansion of the held theory around 
the vacuum conhguration. Equivalently, it corresponds to an amplitude expansion (in powers 
of the “classical held” <f>±) of the 1PI closed-time-path effective action T[0 + ,0_], which is 
dehned in Eq. (2.11) in Ref. |l[. When A0 q is sufficiently large at 770 , or on a time scale for 
rij: to grow to order unity, the perturbative expansion in A breaks down. 

I 11 many inhationary scenarios, condition ([2.39|) does not hold at ijq. A correct analysis 
of the dynamics of the inhaton held must, therefore, be nonperturbative, if the inhaton 
is self-interacting and/or coupled to Bose helds. Again, of interest in “preheating” is the 
time scale for damping of the mean held 0 due to back reaction from particle production 
into the inhomogeneous modes of the fluctuation held. This quantum particle production is 
known to occur by parametric amplification of quantum vacuum huctuations, for the zero- 
temperature, unbroken symmetry system under study here. Boyanovsky et al. [|101|| have 


obtained an approximate analytic expression (in Minkowski space) for the time scale 77 for 
the variance (<^|) to grow to the point where A(</4)/2 is of the same order of magnitude as 
the tree-level effective mass m 2 + A0 2 /6, 


A = 


( (1 + /7 2) 

5(/r l A/B(7)/(87r 2 ) 


-In 


(2.41) 


The function B(f) is of order unity, and in terms of the asymptotic value of / at t] —> 00 , 
B[f(rj —*• cx))] ~ 0.285 953. Their result is valid in hat space and based on a solution of 
the one-loop dynamics which neglects the back reaction of particle production on the mode 
functions. The essential feature of the time scale 77 is that it depends on the ln(A _1 ). As 
a consequence of the analytic solution to the classical mean-held equation and the estimate 
for 77 , it is possible to estimate (for the case of Minkowski space) the effective equation of 
state 7 C for the mean held ||101|| , 


7c = 


Pc 

Pc 


-|/o [ l - i/o 2 ] + |(1 + fo ) 


1 - 


E(k) 

K(k) 


i/ 0 2 1 1 + i/ 0 2 ] 


(2.42) 


where E{k ) is the complete elliptic integral of the second kind ||138|| , p c is the pressure of the 
mean held, and p c is the energy density of the mean held, dehned in Eq. (|2.29| ). The late¬ 
time effective equation of state can be studied using an idealized two-huid model consisting 
of the classical mean-held oscillations and a relativistic component corresponding to the 
energy density in the quantum modes p Q [dehned in Eq. (|3.14|) below]. 

The physical processes discussed above neglect collisional scattering of excitations of 
the inhomogeneous modes due to the A$ 4 self-interaction, for example, binary scattering. 
These scattering processes ultimately lead to thermalization of the system. A quantitative 
understanding of the time scales for such processes in the nonperturbative regime studied 
here within a rigorous held-theoretic framework is at present lacking. A perturbative treat¬ 
ment of collisional thermalization of the system using the Boltzmann equation assumes a 
separation of time scales for collisionless processes ( 77 ) and thermalization. However, due to 
the nonperturbatively large occupation numbers which arise in the resonance band of the 
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inhomogeneous field modes on the time scale 7i, such a naive approach would predict that 
the time scale for thermalization is on the order of (or earlier than) the preheating time scale 
t i. A nonperturbative approach to studying the collisional thermalization of the system is, 
therefore, required. However, within the 1/A expansion (to be discussed in Sec. [fll A| ) , the 
collisional scattering processes are subleading order in 1/A, and thus the separation of time 
scales is assured within this controlled expansion ||101|1 . Let us denote the time scale for 
scattering by r 2 . 


11 


m 


In typical inflationary scenarios, the self-coupling A of the inflaton is very weak 
the range 1CC 12 -10~ 14 (see, however, |55|,|SB[). In our numerical work, / is initially unity, in 
which case the three time scales r 0 , T\ and r 2 separate dramatically, 


Ti/Tq ~ O 


In 


A 


T 2 / 7-1 ~ 0(N). 


(2.43a) 

(2.43b) 


The period leading up to iq is called preheating , because (i) the energy transfer from the mean 
field is entirely nonequilibrium in origin, and (ii) the occupation numbers of the fluctuation 
field are extremely nonthermal. In this regime, since r 2 7q, collisional effects can be 
neglected. In a collisionless approximation, the damping of the mean field is due to energy 
transfer into the inhomogeneous quantum modes, a process similar to Landau damping in 
plasma physics ||119| . 

So far in this section we have not included the effect of spacetime dynamics on the 
particle production and back reaction processes. Cosmic expansion introduces an additional 
time scale // ', where H is the Hubble parameter defined in Eq. (|2.35|) . In typical chaotic 
inflation scenarios, the initial inflaton amplitude can be as large as M P / 3, leading to 


FT 1 ~ 


3t 0 

V2tt 


(2.44) 


at the onset of reheating. In this case, // -C T\ when A is very small. Clearly, for 
sufficiently large initial inflaton amplitude, it is necessary to include the effect of spacetime 
dynamics in a systematic study of preheating dynamics of the inflaton field. 


III. O(A) INFLATON DYNAMICS IN FRW SPACETIME 

In this section, we study the nonequilibrium dynamics of a quartically self-interacting, 
minimally coupled, 0(A) field theory (with unbroken symmetry) in spatially flat FRW space- 
time. We use the covariant evolution equations derived in |l|, in order to study the dynamics 
of the mean field, variance, and the spacetime, at leading order in the 1/A expansion. 


A. The O(A) model in the 1/A expansion 

The classical action for the unbroken symmetry 0(A) model in a general curved space- 
time is 
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(3.1) 


.^ (oY'i,J = -i f d l x^<j 

Z J M 

where the O (N) inner product is dehned byf] 


0 • (□ + m 2 + £R)(j> + ^(0 • 0) 2 


0 • 0 = 0*0 J 0i7- 


(3.2) 


As in Eq. (|2.1| ), A is a coupling constant with dimensions of 1/h, and £ is the dimensionless 
coupling to gravity (and is necessary in order for the quantized theory to be renormalizable). 

In [JT[, the covariant mean-held equation, gap equation, and geometrodynamical held 
equation were computed for this model at leading order in the 1/N expansion. The evolution 
equations follow from Eqs. ( |2.7a|) - (|2.7c|) , with the 2PI, CTP effective action truncated at 
leading order in the 1/N expansion. At leading order in the 1/IV expansion, we need only 
keep track of one component of the CTP two-point function G a b(x, x')\ we choose G ++ (a;, x'), 
which is the Green function with Feynman boundary conditions. The covariant gap equation 
for G ++ at leading order in the 1/N expansion is 

D x + m 2 + £R(x) + ^4> 2 (x) + ^G ++ (x, x) \ G ++ (x, x') = 5 A (x - x')—=L=, (3.3) 

22 J V-9 


plus terms of 0(1/N). The covariant 5 function is dehned in Ref. [|T8l . The mean-held 
equation is, at leading order in 1/N, 


□ + m 2 + (R + lit? + 1 , ix, x) ) <p(x) = 0, 


(3.4) 


where we note that G ++ (x,x) = G a b(x,x ) for all a, b 6 {+,-}■ The coincidence limit 
G ++ (x,x) is divergent in four spacetime dimensions, and the regularization method is de¬ 
scribed in Sec. ms below. The geometrodynamical held equation is 

G, u + A g, u + c + b = -8v tG(T, u ) (3.5) 

in terms of the (unrenormalized) energy-momentum tensor computed at leading order in the 
1/N expansion, which is shown in Eqs. (5.37) and (5.38) in Ref. (lj. 


B. Restriction to FRW spacetime 

Let us now specialize to the spatially hat FRW universe, with initial conditions ap¬ 
propriate to post-inflation dynamics of the inhaton held. As discussed in Sec. m initial 
Cauchy data for </>, G ++ , and a are specihed on a spacelike hypersurface E r;o (at conformal 
time rjo). The spatial symmetries of </> and G ++ for a quantum state \<fi) consistent with a 


5 In our index notation, the latin letters i,j,k,l,m,n are used to designate O (IV) indices (with 
index set {1,..., N}), while the latin letters a, b, c , d, e, f are used below to designate CTP indices 
(with index set {+,—}). 
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spatially homogeneous and isotropic cosmology are given in (|2.21a| -b). As a consequence of 
these symmetries, both the mean held cj) and variance (<^ 2 ) are spatially homogeneous, i.e., 
functions of conformal time only. 

Eq. (|3.3|) for G ++ in spatially hat FRW spacetime has the formal solution 


G++{x,x') = a{rj) 


j \-1 


d 3 k 

( 2 nj* 


D ik-(x-x') [Q 


[©(V - v)uk{v)*u k W) 


+©(?/ - r)')u k {r)')u k (rj)% (3.6) 


in terms of conformal-mode functions u k which satisfy a harmonic oscillator equation with 
conformal-time-dependent effective frequency, 


(A + % = °. 


(3.7) 


The fact that Ukiv) depends only on )] and k (where k is comoving momentum) implies that 
G ++ is invariant under simultaneous spatial translations and rotations of x and x!. The 
effective frequency Q k {v) appearing in Eq. O is defined by 


ni( V )=k 2 + a 2 m\r ] ), 

(3.8a) 

m 2 (v)=m 2 (v)+(z-^) R(v), 

(3.8b) 

M 2 ( V ) = m 2 + ^ 2 (v) + 

(3.8c) 


Initial conditions for the positive frequency conformal mode functions u k {ji) must be specihed 
(for all k) at r]o. A choice of initial conditions corresponds to a choice of quantum state 10) 
for the fluctuation held <p H ; initial conditions are discussed in Sec. |1V A| below. The (bare) 
variance (<p|) has a simple representation in terms of the conformal-mode functions: 


{<,M<Pn(x) 2 \4 >} = hG ++ (x,x ) = 


h 


d 3 k 

(2 


\Mv)\' 


(3,9) 


It should be noted that this expression is divergent, in consequence of our having computed 
the variance in terms of the bare (unrenormalized) constants of the theory. In terms of 
a physical upper momentum cutoff K, G ++ (x,x) diverges like A' 2 ; there is additionally a 
logarithmic dependence on K. In addition, the mode functions Uk depend on (ip |) through 
Eq. (|3.8c|) . The leading-order, large- N, mean-held equation in spatially hat FRW spacetime 
becomes 

„ On' „ 

(j)" H- (j)' + a 2 M 2 (rj)(j) = 0, (3.10) 

a 

where the time-dependent bare effective mass M{rj) is given by Eq. (|3.8cj) . For simplicity of 
notation, we will henceforth write M instead of and similarly for 9R(r)). 

Finally, we can express the bare energy-momentum tensor in terms of the conformal-mode 
functions u k {rj). As discussed in Sec. E3 it is convenient to work with the 00 component 
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and the trace of the energy-momentum tensor. The components of the classical part of the 
energy-momentum tensor are spatially homogeneous, and given by 


T 0 C 0 (V ) 

r c (v) 




2 


2 


+ 


2 a' 


+ 1 a 2 (m J + ^ 


+ 


3£(a 


f \2 


1 f (65-i)(,y) 2 +6«0+— ^ 


+ 2 m + 


2a 4 
A , 


> 2 + -i? 


(3.11a) 

(3.11b) 


The quantum energy-momentum tensor components are also spatially homogeneous. We 
find for the 00 component, 


TS,M = 


h f d 3 k 


2 a 2 J ( 27 t ) 3 


\<\ 2 + (k 2 + a 2 M 2 + (l-60^)|fi t | 


+ (6£ - 1)- {u' k Yu k + u' k ul 
a 


and for the trace, 


? Q (V) = 


h f d 3 k 


(2tt) 


(6e- 


a 


\u' k \ 2 — (k‘ 2 +a 2 M 2 )\uk\‘ 2 — — (■ u' k )*u k + u k u k 
,/\2 


.a 


a l 

,2 D A|~ 12 


a z 


ea 2 l? |h fc | 2 +a 2 M 2 |h fc | 


(3.12) 


(3.13) 


It can be shown by asymptotic analysis that, in terms of a physical upper momentum cutoff 
K , the bare is quartically divergent, i.e., 0(K 4 ), and that (for minimal coupling) T Q is 
quadratically divergent. In addition, the components of the bare energy-momentum tensor 
contain the effective mass M 2 , which contains the divergent variance (</?|). The energy 
density p Q of quantum modes of the tp held is defined in terms of T® 0 by 

P« = 4 (3.14) 

a z 8 

We shall also refer to p Q as the energy density of the “fluctuation held.” 


C. Renormalization of the dynamical equations 


The variance (<^|) and quantum energy-momentum tensor components Tq 0 and T Q are 
divergent in four spacetime dimensions, and must be regularized within the context of a 
systematic, covariant renormalization procedure. In the “in-out” formulation of quantum 
held theory, renormalization may be carried out via addition of counterterms to the effective 
action, which amounts to renormalization of the constants in the classical action ||140|| . The 


closed-time-path formulation of the effective dynamics is renormalizable provided the theory 
is renormalizable in the “in-out” formulation |32| , |68| , as is the case with the O (N) held theory 
in curved spacetime [jiyi41|,|T4^j. For our purposes it is convenient (in this model) to carry 


out renormalization in the leading-order, large-iV, evolution equations, rather than in the 
CTP effective action |67I|. 
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In this study we employ the adiabatic regularization method of Parker, Fulling, and Hu 
The idea is to define an adiabatic approximation to the conformal mode function, 


102,103 


and then to construct a regulator for the integrands of the bare energy-momentum tensor 
and variance from the adiabatic mode functions EH' Renormalization occurs when we 
define the renormalized variance and energy-momentum tensor to be the difference between 
the bare expressions and the regulators and simultaneously replace the bare quantities m, 
A, G, b , c, A, and £ by their renormalized counterparts. The equivalence of this procedure 
to other manifestly covariant methods (such as dimensional continuation) is well established 
143|| . We implement renormalization as a two-step process: First, we adiabatically regularize 
the variance and renormalize £, m, and A; Second, we adiabatically regularize the energy- 
momentum tensor and renormalize the semiclassical geometrodynamical held equation. 

We define the adiabatic order of a conformal mode function as follows: let 12^ (77) 
h2fc(r//T), where T is introduced as a time scale which is formally taken to be unity at the 
end of the calculation. Then the adiabatic order of an expression involving derivatives of 
Hfc is simply the inverse power of T, of the leading-order term in an asymptotic expansion 
about T —> oo. However, in order for the adiabatically regulated energy-momentum ten¬ 
sor for an interacting scalar held theory to agree with the renormalized energy-momentum 
tensor obtained by manifestly covariant methods (e.g., covariant point splitting ||131|| ), it is 
necessary to define the adiabatic order of expressions involving A and derivatives with re¬ 
spect to 7 /, such as A(0 2 )", as the sum of the exponent of <j) and the number of conformal time 
differentiations ||144|| . Therefore, A((^|) // is considered fourth adiabatic order, as is A(0 2 )". 

Having dehned adiabatic order, we now construct the adiabatic mode functions. It is 
well known that the WKB ansatz 


Mb) = 


V'iwW) 


exp i i 


dr\W (rj'] 


(3.15) 


turns the harmonic oscillator equation (with time-dependent frequency 12*,) into a nonlinear 
differential equation for W, 


Mb ) 2 = n 2 k ( v ) + 


3\W'(rj)] 2 W"(rj) 
4 W 2 (r]) ~ 2HM)' 


(3.16) 


Starting with the lowest-order ansatz W^\rf) = f4(b)> one can iterate this equation; the 
nth-order iteration yields the nth-order WKB approximation for Uk . For the free held theory, 
the nth-order WKB approximation gives an expression for Uk which is of adiabatic order 2n. 
In the interacting case, the above definition of adiabatic order calls for removing terms such 
as \((j) 2 )"" at 4th adiabatic order. Thus we have a method of deriving expressions for T® 0 , T q , 
and (y?!) at fourth, fourth, and second adiabatic orders, respectively. One then sets T = 1 in 
the truncated expression. We can thus obtain a fourth-order adiabatic approximation to the 
quantum energy-momentum tensor (T^) ad4 , and a second-order adiabatic approximation 
to the variance (</? H ) ad2 . By subtracting (T® u ) ad4 from the divergent T® u and (<^ 2 ) a d 2 from 
the divergent (<^|), finite expressions for the renormalized energy-momentum tensor and 
variance are obtained. 

First we regularize the variance (< p 2 ), and carry out a renormalization of A, m, and £. In 
the leading-order, large -N approximation, no terms appear in the mode-function equation 
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( p.7| ) which would permit addition of counterterms; therefore, Q k must be finite ||145|| . The 
effective frequency Q k which appears in Eq. (|3.8a ) is the “bare” effective frequency, which 
we denote by (fl k ) B - In conjunction with the adiabatic regularization procedure, we fix the 
renormalization scheme by demanding equivalence of the bare and renormalized effective 


mass 


119 


Mh = (^) 




(3.17) 


where the “R” subscripted quantities are renormalized, 
have 


Using Eqs. (|3.8a|) and ( |3.8b|) , we 


£rR + — £ b R + M|, (3.18) 

where M| is defined in Eq. ( |3.8c|) , 

Ml = m 2 + y0 2 + y(v?h)b, (3.19) 

and Ml is defined similarly, 

M b = m 2 + —ft + —{vV)*.- (3.20) 

Now, A b , m B , and £ B are the bare constants of the theory which appeared (without B’s) in 
the classical action ( |3.1| ). The renormalized quantities in Eq. (|3.8a|) are dehned below. The 
bare (<^|) B is a conformal-time-dependent function dehned by Eq. ( |3.9|) , 

(vI(v))b = y/ (^3^ fe ^| 2 ’ ( 3 - 21 ) 

where the conformal-mode functions u k {rf) obey Eq. (0- Now we demand that the renor¬ 
malized theory be minimally coupled, i.e., we set £ R = 0. Because of Eq. ( |3.17|) , we can 
formally use (0^) R in computing the adiabatic regulator for the variance (</j|) B . Computing 
the asymptotic series (in 1/T) of the quantity \u k (r])\ 2 to 0(1/T 2 ), where h2 2 (?//T) is the 
effective frequency, we obtain (after setting T = 1) 



ZO 


d 3 k M (C') 2 — 2CC'" MIC" 5M^(C') 2 ' 

(27t) 3 _cu k 8C 2 ivl + 8u>l 32 u 7 k 


in terms of auxiliary functions 


C(v) = a 2 (^), 


and 


D{n) 


C'(y ) 

c(v )' 


In Eq. (|3.22|) the symbol t o k is dehned as follows 


(3.22) 


(3.23) 


(3.24) 
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ul = k 2 + a 2 Ml. (3.25) 

In the adiabatic prescription, the renormalized variance (</?|} R appearing in Eq. ( |3.20| ) is 
defined by 

(^)h = (^h)b - <v!U (3-26) 


where the first term on the right-hand side is given by Eq. ( |3.21| ), and the second term on the 
right-hand side is given by Eq. ( [3.22| ). Everything on the right hand side can be expressed 
in terms of renormalized quantities, so this procedure is well defined. Written out explicitly, 
the renormalized variance satisfies the equation 


(A = 


h 


d 3 k 


C J (2vr) s 


\u k \ 2 ~ 




(C") 2 - 2CC" 
16 C 2 Col 


+ 


MIC" 

16 uj\ 


5M*(C") S 
64 Col 


(3.27) 


One can use the WKB approximation for u k {r /) to compute the asymptotic series for the 
integrand in Eq. (|3.27|) in the limit k —> oo, and show that the integral is convergent. Since 


M'l is contained in the integrand above, Eq. (|3.27|) leads to an integral equation for the 


renormalized effective mass M, 


R> 


7 i r2 2 A R ?2 ^lA R f d k 

Ml = ml + -^(b 2 + — I - 


2 C J ( 2 tt ) ; 


\u k \ 2 - 


2 uii 


{C') 2 - 2 CC" MIC" 


5Mi(cy 


1 6C 2 Co 3 


16o >1 


64 


(3.28) 


Eqs. (|3.27|) and (|3.18|) together define A R and m R . All physical quantities should now be ex¬ 
pressed in terms of the renormalized parameters m R and A R of the theory. The renormalized 
mean-field equation now reads 


</>" + —0' + a 2 M 2 0 = 0, 

a 

where M 2 is given by Eq. (|3.28| ), and the mode functions in Eq. 
equation, 


(3.29) 

obey the homogeneous 


w- +k2+a2Ml 


u k = 0. 


(3.30) 


rvA 


The initial conditions for the conformal-mode functions at r/o are discussed in Sec 
below. 

Having obtained a renormalized mean-field equation, we now turn our attention to reg¬ 


ularizing the quantum energy-momentum tensor. As a consequence of Eq. (|3.17|) , we can 


substitute M —> M R and £ —s► £ R in the equations for the components of the quantum energy- 


momentum tensor, Eqs. (|3 .1 % |3.13|) . Since we wish to study the minimal coupling case, we 
set i^ R = 0. To avoid confusion we denote the bare energy-momentum tensor components 
( p.l2| ) and (|3.13|) by (T 0 ^) B and (T Q ) B , respectively. Let us also relabel the bare constants 
b , c, G, and A appearing in the bare semi classical Einstein equation (|3A| ) as b B , c B , G B , and 
A b . Applying the method described above to construct the adiabatic regulator, for Tq 0 we 
find 
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PXL« = 


ft 


d 3 k 


AC J (2t r) 3 


(cy 

2uJk+ 4 C++ 

r2\l 


Ml(C') 2 9 (C") 


/\4 


AC 

i\in" 


c\Miy (< eye" , (c ") 2 c'c 


+ 


+ 


64C' 4 


w 


4 ' 4C 3 16C 2 8C 2 

^1?" + 3^3(“ 5(C "1 4 + 4C\C')(Ml)' 


+io cyc'fc" + 2 c\c"f - ac 2 c'c"' 


J a;; 4 


YALi 


-5(C") 4 + 40(C") 2 C"' + 2C 2 (C") 2 
1 7Ml(C') 2 / __ _\ 1 


V128C' 2 / 

-AC 2 C'C"'] i + 7iV Jg^ 2 (-5(C") 2 + 2C'C" / ) 


vi 




105M|(C") 4 1 


1024 


v 


li r > 


where Vk is defined in Eq. fl3.25|) . For T, we find 


(T Q ) ad4 = 


h 


d 3 k 


2 C 2 J (2t r) 3 


M 2 C-^! + ^U 
R 2C 2 + 2(7/ Uk 

M 2 / \ 1 

L8c(“ 3(C,,)2 + 4CC,/, ) + - 


9(C') 4 + 4C 4 C"(M, 


16C7 4 

-21C'(C , ) 2 C" / + 6 C 2 (C") 2 + 4C' 5 (M 2 )" + 8 C 2 C'C'" 


t n/r2y 


-2 C 6 C 


3 rim' 


VI 


Ml 
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/\4 


Mi 


(-3 (C") 2 + CC") 


- 64C' 4 C"(M 2 ) / - 208C'(C") 2 C 


i\2riii 


128 C 3 

+60 C 2 (C") 2 + 16C 5 (M 2 ) // + 80C' 2 C"C"" - 16C' 3 C"" / 


J vi 


+ 


5CMI 


Mf 

32 v ~ ' ' 3’2C 2 
40C'(C") 2 C" / + 15(C") 2 (C") 2 + 20 C 2 C'C'" - C 3 C"" 


“(C+ + 3+5(l5(C+ - 10C 4 C'(M 2 )' 


J Vi 


(3.31) 


+ 


7 Ml 
L256 C 
21 M|(C") 2 


(l5(C") 4 - 80C(C') 2 C" + 6C' 2 (C") 2 + 8 C 2 C'C"'^ 

(15(C') 2 -11CC'')4 + 1155CM “ ' 

v / vi 


vi 


2048 


(CT+ia 


(3.32) 


In the free-field limit (A R = 0), the regulators (|3.31|) and (|3.32|) agree with the minimal- 


coupling, spatially flat limit of the adiabatic regulators obtained by Bunch 143 


The renormalization procedure for the semiclassical Einstein equation ( |3.5| ) can now be 
precisely stated. According to the adiabatic prescription, we define the quantum energy- 
momentum tensor by 


(+„)» = CC)b - (7?„W 


' fns 


ins ) 


(3.33) 
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It can be checked that the momentum-integral expressions for the two independent compo¬ 
nents of Eq. (|3.33) are convergent. In terms of the total energy-momentum tensor 

(after renormalization) is 

(T„„), = (T£)» + (T« ) R - A ((rf)» ) V. (3.34) 

where (T^) R stands for T° v , and renormalized quantities are substituted for bare quantities. 
The bare quantities G B , A B , b B , and c B are now replaced by G R , A R , b R , and c R in the 
renormalized semiclassical geometrodynamical held equation, 

+ Aj+ c R + b R ^ H^ v = —8(/r G R (r ixv ) R . (3.35) 


D. Renormalized semiclassical Einstein equation 


Using semiclassical methods to study the dynamics of the inhaton held in FRW spacetime 
requires that the Hubble parameter be much less than the Planck mass, H <§C M P . On 
dimensional grounds, c R and b R are likely to be of order h 2 M R ~ 2 , in which case R c r R 2 , 
and R b R R a/3 R a p, provided R a p ^ 0. Let us, therefore, set b R = 0 and c R = 0, and 
additionally, let us choose A R = 0, so that Eq. (|3.35|) becomes the renormalized semiclassical 
Einstein equation (without cosmological constant), 


G^ v — SitGk 


CC). + CC)»-^(<rf>«) 


(3.36) 


Taking the trace of Eq. ( |3.36|) in spatially hat FRW spacetime, we find 

Ah 


% = 8ttG r 
a 6 


(T c ) r + (T«) r -^((^) r ) 2 


(3.37) 


Recalling that £ R = 0, and using Eq. ( |3.11b|) , the classical part of the trace of the renormal¬ 
ized energy-momentum tensor is given by 


(^°)k = 


-(0O 2 + 2 ( ml + A 4 


R 12 


(3.38) 


and the quantum trace of the renormalized energy-momentum tensor is given by 
d 3 k 




(2tt) 3 _ 


\u 


/ 2 


(k 2 - 2a 2 M 2 )|« fe | 2 - - [«)*«„ + +LA|fi t | 2 


(T Q )^, (3.39) 


where (T Q ) ad4 is dehned in Eq. (|3.32|) . As discussed in Sec. m the 00 component of the 
semiclassical Einstein equation is a constraint, and is given by 


3 (a' 


/\2 


= 87 tG t 


(IS). + (IS). - 


(3.40) 
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From Eq. (|3.11a|) , the expression for the classical part of the 00 component of the renormal¬ 
ized energy-momentum tensor is given by 


1 ' V\2 , 1 2 


(T 0 C o)r = x(0') J + 


2 , ^R 12 

m H + -J<t> 


(3.41) 


and the quantum part of the 00 component of the renormalized energy-momentum tensor 
is given by 


(^oo)r — 


h r d 3 k 

2 a? J ( 2 vr ) 3 


Wk\ 2 + (^ 2 + a 2 M 2 ) \u k \ 2 — — [(u' k )*Uk + u' k u 


(^Oo)ad4, (3.42) 


where (T^)^ is defined in Eq. (|3.31|) . 

Eqs. (13.37|) and (|3.29|) are coupled differential equations for a and 0, involving complex 
homogeneous conformal-mode functions u k which satisfy Eq. ( j3.30| ). The conformal mode 
functions are related to the variance (<p 2 ) R by Eq. (|3.27|) . This is a closed, time-reversal- 
invariant system of equations. The initial data at r] 0 must satisfy the constraint equation 
( |3.40| ). We now drop all “R” subscripts, because we will henceforth work only with renor¬ 
malized quantities. 


E. Reduction of derivative orders 


The adiabatic regulators (|3.22|) , ( |3.31|) , (|3.32|) for the variance and energy-momentum 
tensor contain derivatives of up to fourth order in a and up to second order in 0 2 and (< p\ ). 
The presence of the former can be understood as resulting in part from the well-known trace 
anomaly for a quantum held in curved spacetime [ 146 1, which contains higher-derivative 
local geometric terms, e.g., OR. In addition, there are nonanomalous finite terms which 
result from the renormalization of the energy-momentum tensor and the choice of minimal 
coupling. 

The effect of higher derivatives in the semiclassical Einstein equation has been much 
studied in the literature ||147| - |151[j . The higher-derivative evolution equations for a and <ft 
have a much larger solution space than the classical Einstein and mean-held equations, and 
in general, the higher-derivative semiclassical Einstein equation is expected to have many 
solutions which are unphysical. In addition, the semiclassical Einstein equation (which is 
fourth order in a) requires more initial data than the classical Einstein equation in order to 
uniquely specify a solution. However, Simon and Parker ||150|jl51|| , following the methods of 
Jaen, Llosa, and Molina [|152|| , have shown that in one-loop semiclassical gravity, there exists 
a procedure for consistently removing the unphysical solutions within the perturbative ( h ) 
expansion in which the equations are derived. The procedure corresponds to the addition 
of perturbative constraints, thereby yielding second-order equations which require the same 
amount of initial data as does the classical Einstein equation. Their method involves re¬ 
ducing the order of the a'” and a'”' terms in the semiclassical Einstein equation using strict 
perturbation theory in h. 

In this study we follow the approach of Simon and Parker to reduce the order of the 
equations for 0, a, and (<p 2 ) to second order. We replace all expressions involving o'" and 
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a'"' with expressions a!" and a"" obtained from the classical Einstein equation, i.e., Eq. fl3.37|) 
with h = 0. This procedure is physically justifiable in this model for the following reason: 
At early times, the dominant contribution to the energy-momentum tensor is classical, T^ u . 
Therefore, the deviations a'" — a'” and a"" — a"", which are entirely quantum in origin and 
oc h, are at early times expected to be very small. In addition, at late times the Universe 
is expected to become asymptotically radiation dominated, in which case a"' = a'"' = 0. 
The classical approximations to the late-time behavior of a and a"" should also have this 
property, regardless of whether the mean-field oscillations are harmonic or elliptic. This 
procedure is, therefore, physically justifiable in the system studied here. 


IV. ANALYSIS 


Having derived coupled dynamical equations (|3.29|) , (|3.37|) , ( |3.30| ) for the mean field <f>, 
scale factor a, and conformal-mode functions Uk, respectively, we now proceed to solve them. 


A. Initial conditions 


At the Cauchy hypersurface at r/ o, we specify initial conditions on the conformal-mode 
functions Uk which correspond to a choice of quantum state for the fluctuation field ^ H - 
Based on the analysis in Sec. 11C . we choose boundary conditions at r/o which correpsond 
to the adiabatic vacuum state for cp H at r] —> —oo. From the semiclassical Einstein equation 
( 12.10| ), the slow-roll condition (|2.33|) , the potential-dominated condition (|2.32|) , and assuming 
that the variance (ipl) satisfies 




(4.1) 


for '// < 7]o, it follows that the spacetime is asymptotically de Sitter at conformal-past infinity. 
Using the approximate solution ( |2.34| ) for the scale factor for r] < r] o, we can solve the mode 
function equation (|3.30|) for rj < r ] 0 at the same (0th) adiabatic order. The general solution 
is of the form 


u k (v) — 


- Tjo) 


c \H0{kCn - H-Hri) - r,„)} + - H~\ r,) - 

(4.2) 


where H B) and H ® are the Hankel functions of first and second kind, respectively ||138| 
and v is defined by 


2 9 

V — - 

4 


M 2 
JP ' 


The function H{rj) is defined as in Eq. (^.35|), 


H(v) = 


8ttGp c 


(4.3) 


(4.4) 
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where now p c = cl 2 Tq 0 . The Hubble parameter must be slowly varying for this approximation 
to hold, i.e., the expansion rate nonadiabaticity parameter [f|^] 

to H = §2 « L ( 4 ‘ 5 ) 

The Wronskian condition on the mode functions (which comes from the canonical held 
commutation relations) requires that 

I4| 2 +I4 | 2 = l- (4.6) 


By choosing c[. and c^, different vacua are obtained. The Oth-order adiabatic vacuum 
(matched at r / = —oo) is constructed by choosing c\ and cf so that Uk smoothly matches the 
positive-frequency Oth-order WKB mode function at 77 = — 00 . This corresponds to c\ = 1 
and c\ = 0, for all k. Using the asymptotic properties of the Hankel function, the adiabatic 
limit k, \t)\ —> 00 , can be derived, and verihed to have the correct form, 


lim Uk — 

A:, 177! —>00 



(4.7) 


In addition, the high-momentum, hat-space limit (k, H^ 1 —> 00 ) gives the same result. The 
initial conditions for the Uk at ij 0 are then defined by demanding that the Uk functions 
smoothly match the approximate adiabatic mode function solutions (for 77 < 770 ) at 7] = tj 0 . 
This leads to the following initial conditions for the conformal-mode functions: 


Mho) = 



u'k(Vo) 


d 

di] 




(4.8a) 

(4.8b) 


where H$ = H( 770 ). The above conditions are valid only at 0th order in the above-defined 
adiabatic approximation, where terms of order H'/H are discarded. It is straightforward to 
show that Eq. ( |4.5| ) is valid given the slow-roll ( |2.33| ) and inflation ( |2.32| ) assumptions. In 
addition to the initial conditions for Uk at 770 , we may freely choose initial values for 0 ( 770 ) 
and < 7 /( 770 ), subject to the constraint that 0' must be small enough that conditions (|2.32| ) 
and (|2.33|) are valid. We are already assuming that < 2 ( 770 ) = 1. The initial value of a\7] 0 ) is 
fixed by the constraint equation (|3.40|). 


B. Numerical solution 


I 11 this section we describe the methods we used to solve the coupled evolution equations 
for 0 [Eq. (|3.29|) 1, a [Eq. ( fi.371 )1, and Uk [Eq. ( p.30|) l numerically.^ We evolved a repre¬ 
sentative sampling of mode functions Uk for the region of momentum space 0 < k < Ka, 


6 Henceforth, we set h = 1 and work in units of energy where m = 1. 
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where K is a physical upper momentum cutoff. [] Employing a physical ||153|| , as opposed to 
comoving, momentum cutoff is necessary because a comoving cutoff would require the use 
of the renormalization group equation to track how the renormalized parameters flow as the 
scale factor a increases at each time step. For a comoving cutoff the quadratic divergence 
in the variance would be proportional to 1 /a 2 , requiring a time-dependent renormalization 
(see mi, for example). The use of a physical upper momentum cutoff yields a quadratic 
divergence which can be removed by a non-time-dependent mass renormalization [ |153| . 

We chose a variety of values of K/m between 50 and 70. The sampling of momentum- 
space is carried out with a uniform binning, with total number of bins N bina . Various values of 
N bina were used, all greater than 10 4 . Eq. ( |3.28|) was solved by iteration, and the momentum 
space integrations were performed numerically using the 0(1/N bins ) extended Simpson rule. 
The differential equations ( p.37| ) and ( |3.29| ) were evolved using 4th-order Runge-Kutta with 
adaptive step-size control; the target precision for the time steps varied between 10 -6 and 
10” 8 . Cutoff independence was verified a posteriori by explicitly checking that the results of 
the numerical solution were insensitive to a rescaling of K/m. The solutions were computed 
to a conformal-time scale of 400 m" 1 . A typical solution computed according to the above 
methods required on the order of 300 h of CPU time on a modern workstation. 


C. Results 


A primary goal of this work is the quantitative study of the effect of spacetime dynam¬ 
ics on the parametric resonance energy-transfer mechanism in nonequilibrium zero-mode 


oscillations of a quantum held. As discussed in Sec. [II Dj . this energy transfer, and the corre¬ 
sponding damping of the mean held due to back reaction, occur on a time scale of order T\ 
dehned in Eq. (|2.41|) . We numerically evolved the evolution equations for a, (p, and (<p 2 ) for 
various values of M P /m, ranging from very large values (corresponding to Minkowski space), 
to small values (corresponding to a strong-curvature, rapid-expansion regime). Figs. 0-0 
show the resulting time dependences for the mean held (p , the scale factor a, the variance, 
A(<p 2 )/2, the energy density p, the energy density in quantum modes p Q [dehned in Eq. 
Q3.14Q ], and the pressure-to-energy-density ratio 7 . The different solutions plotted corre¬ 
spond to different values of M P /m, with A = 10” 14 , K/m = 50, and cp(r] 0 )/m = 2.0 x 10'. 
As discussed in Sec. |IV B[ a physical momentum cutoff K was used. The values chosen for 
Mp/m were 10 14 , 10 12 , 6 x 10 10 , and 6 x 10 9 . The choice of 0(p o ) and A hxes r/o by Eq. 
( |2.34| ) and H 0 by Eq. ([2.35| ). Table | shows the values of M P /m, the inverse Hubble constant 
i7 - 1 (p 0 ), and the hgure numbers in which the corresponding solutions are plotted. 

The time scales dehned in Sec. [II 1)| can now be explicitly computed. Using Eqs. (|2.38a|) 


and (|2.35|) , we have /( 70 ) = \/2, Po = 0o, and 


H(V 0 ) = 


87 T(p\ 


0 


3M P 2 


= Ho 


(4.9) 


' A finite momentum cutoff is necessary due to the triviality of the theory when the cutoff is taken 
to infinity [119]. 
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Using Eq. (|2.37|) , we find r 0 ~ 4.118 32 m~ l . The value of T\ predicted by Eq. (|2.41|) is 
132.624 m _1 , which is very close to the value predicted by Eq. ( 2.43a ), 132.759 m~ l . For 
the cases M P /m = 10 14 and 10 12 , it is clear from Table | that Hq 1 r 1; so that the effect of 
cosmic expansion is expected to be insignificant on the preheating time scale 17 . For the case 
M P = 6 x 10 10 , 1/(H 0 ti) ~ 7.8, so that the effect of cosmic expansion should be apparent 
and non-negligible. For the case M P = 6 x 10 9 , 1/(HqT\) ~ 0.78, and cosmic expansion 
should have a significant effect on parametric amplification of quantum fluctuations. 

Figs. M show the dynamics of the mean field and variance in the regime of very weak 
cosmic expansion, H~ l C q. As expected, under the influence of the clliptically oscillating 
mean field, the variance (<£> 2 ) grows exponentially in time until A(</ 3 2 )/2 is of the same order 
as m 2 + A0 2 /2, at which point back reaction shuts off the resonant transfer of energy to 
the inhomogeneous modes. The time scale for the variance to become of order unity can be 
clearly seen to be ~ T\. As seen previously in studies of preheating dynamics in Minkowski 
space [TUI], on the time scale ~ r 1; the mean field decouples from its own fluctuations and 
oscillates with an asymptotically finite amplitude, given by ||101|| A0 2 /(2m 2 ) = 0.914. In the 
Minkowski space limit, corresponding to M P /m —> oo, covariant conservation of the energy- 
momentum tensor implies that dp/dtt = 0. This was verified for the case of M P /m = 10 14 , 
where no change in p was detected to within the numerical precision of our algorithm, as 
expected from dimensional analysis of Eq. ( |3.37| ). The increase in the scale factor for these 
cases was within a few parts in 10 6 of the initial value a(^ 0 ) = 1. The asymptotic equation of 
state plotted in Fig [13| is observed to be 7 ~ 0.18. This is exactly what would be predicted 
for a two-fluid model consisting of a mean field with equation of state given by Eq. ( [2.42|) , 
7 C ~ 0.0288, and a relativistic gas corresponding to the energy density of the ip held, with 
7 q ~ 0.3333. The average 7 Q + 7 c = 0.182. 

For the case M P /m = 6 x 10 10 , the effect of cosmic expansion is clearly visible in Figs, p- 
T4| . In Fig. |9j. the coherent oscillations of the mean held for the time period 0 < r] — r/o < 
~ 27 t 0 are clearly seen to be redshifted by the usual 1/a factor expected from the Hubble 
damping term in Eq. ( |3.29| ) . The expected asymptotic equation of state (taking into account 
cosmic expansion) computed from a simple two-huid model is ~ 0.133, in agreement with 
Fig. [H| 

Figs. show the solution for M P /m = 6 x 10 9 . In this case, 1/(H 0 ti) ~ 0.781. From 

Fig. [L7|, we clearly see that cosmic expansion renders parametric amplification of quantum 
huctuations an inefficient mechanism of energy transfer to the inhomogeneous modes. The 
very rapid oscillations of the mean held at late times are due to the conformal time scale used 
here, in which the oscillation period of the mean held decreases inversely with a. Damping 
of the mean held due to cosmic expansion is the dominant effect in Fig. [15]. The power- 


law decrease in energy density consistent with matter having an effective equation of state 
7 ~ 0.0288 can be seen in Fig. [18]. At 7 = 300 m -1 , the ratio p Q /p ~ 0 . 0002 , so the fraction 
of energy density in the inhomogeneous modes is negligible in comparison to the classical, 
mean-held contribution. Since the variance (<^|) is never large enough that it dominates the 
effective mass M, the mode functions approximately obey the one-loop equation, in which 
the effective frequency is k 2 + a 2 {m 2 + A0 2 ), neglecting the a”/a term. The width of the 
resonance can then be shown to be approximately given by k 2 < A</>q/ 2. The variance is 
damped by 1/a 2 due to cosmic redshift, so when H~ 1 ~ 77, the variance never grows to be 
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of order unity. 

In addition to varying M P , the coupling A was varied, with results in agreement with Eq. 
(p.41|), showing a logarithmic dependence of T\ on A -1 . 


V. DISCUSSION 


The main objective of this series of papers is to investigate the nonperturbative and 
nonequilibrium quantum processes in the reheating epoch of inflationary cosmology. In this 
paper we use the minimally coupled, quartically self-interacting scalar O (N) field theory 
as a model for the inflaton field, and study its nonequilibrium dynamics nonperturbatively 
in a spatially flat FRW spacetime whose evolution is driven by the quantum field. We 
solve the coupled, self-consistent semiclassical Einstein equation, mean-field equation, and 
conformal-mode-function equations numerically. Our goal in this paper is to study the effects 
of spacetime dynamics on the mean field, and parametric amplification of quantum fluctu¬ 
ations. This process of energy transfer from the mean field to the inhomogeneous modes is 
inherently nonperturbative and nonequilibrium. It requires the use of the closed-time-path 
formalism and the two-particle-irreducible effective action. As our focus in this paper is on 
the parametric amplification of quantum fluctuations, we assume unbroken symmetry. Our 
analysis is, therefore, most relevant to reheating in chaotic inflation scenarios. We use the 
two-loop, covariant equations for the mean field and the two-point function for the fluctua¬ 
tion held derived in the preceding paper []1[ , and study the case of leading order in the 1/N 
expansion, an approximation which is valid on time scales much shorter than the mean-free 
time for multiparticle scattering ( 72 ). For FRW spacetimes, we use the well-established 
adiabatic regularization procedure to obtain finite expressions for the renormalized variance 
and energy-momentum tensor which enter into the mean-held equation, conformal-mode 
function equations, and the semiclassical Einstein equation. In our approach, covariant 
conservation of the energy-momentum tensor is preserved at all times, as it should be. (It 
should not and need not be put in by hand, as was done in a recent study of reheating 
in a hxed background FRW spacetime ||129|| .) We use the adiabatic vacuum construction 


(matched at conformal-past infinity) to dehne the quantum state for the fluctuation held 
in FRW spacetime with asymptotic de Sitter initial conditions; this is the most physical 
vacuum construction given the decidedly nonadiabatic conditions which prevail at the end 
of inflation. The instantaneous Hamiltonian diagonalization constructions used in earlier 
studies of reheating in curved spacetime [|12(3|J118|| are known to be problematic ||133| . 


We evolved the coupled dynamical equations for the mean held, variance, and scale factor 
using standard numerical methods, for time scales of 400 m -1 , where the initial period of 
mean-held oscillations is 4.11832 m -1 . Several regimes for the parameters of the system 
were investigated. From the solutions of the dynamical equations we studied the behavior 
of the scale factor a, the mean held 0, the energy density p, pressure-to-energy-density ratio 
7 , and the inhomogeneous-mode (fluctuation-held) energy density p Q . The solutions of the 
dynamical equations were analyzed for a variety of values for M P ro, the parameter which 
controls the rate of cosmic expansion relative to the time scale for mean-held oscillations in 
the model. In the case of negligible cosmic expansion, corresponding to very small initial 
inhaton amplitude, the dynamics is identical to that seen in the group 2B (see Sec. jPA[ ) 
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studies of 0 (N) preheating in Minkowski space [|101|| . In particular, the conservation of 


energy and logarithmic dependence of the preheating time scale iq on the inverse coupling 
A -1 [as shown in Eq. ( |2.41|) ] are confirmed. For the case of moderate cosmic expansion, 
H(rj 0 )ti 10 [where H(r ) 0 ) is the Hubble parameter at the initial time ?7 0 ]> energy transfer 
via parametric amplification of quantum fluctuations is still efficient, and the dynamics can 
be understood using the analytic results of [|10 1|| (for Minkowski space), in terms of the 

conformally transformed mean-field amplitude (j) = (j)/a and oscillation period fo = To/a. 
The asymptotic effective equation of state is found to be consistent with the prediction of a 
simple two-fluid description of the late-time behavior of the system. 

The most significant physical result concerns the case of rapid cosmic expansion, where 
H~\r, 0 ) — T\. In this case we find that parametric amplification of quantum fluctuations 
(via parametric resonance) is an inefficient mechanism of energy transfer to the inhomoge¬ 
neous modes of the inflaton, because the parametric resonance effect is inhibited both by 
redshifting of the mean-field amplitude and by the redshifting of the physical momenta of 
the modes out of the resonance band. The energy density of particles produced through 
parametric resonance is in this case redshifted so rapidly that, in our model, the term 
A(<^!)/2 never grows to be of the order of the tree-level effective mass, m 2 + A</> 2 /2. As 
the mean-field amplitude is damped (oc 1/a) due to cosmic expansion, eventually the res¬ 
onant particle production ceases, and the mean field oscillates with a damped envelope at 
late times. This leads us to the following conclusions: (i) On the physical level, in chaotic 
inflation scenarios with a A<h 4 inflaton minimally coupled to gravity and with a large initial 
inflaton amplitude at the end of slow roll, parametric amplification of the inflaton’s own 
quantum fluctuations is not a viable mechanism for reheating the Universe, unless the self¬ 
coupling is significantly increased.^ This does not imply that the phenomenon of parametric 
amplihcation of quantum fluctuations does not play a vital role in the “preheating” period 
of inflationary cosmology, for different models and/or couplings. The interesting case of a 
<f> 2 X 2 model will be addressed in a future publication |$6|. (ii) On a more methodological 


level, we conclude that a correct study of the reheating period in a chaotic inflation model 
with large inflaton amplitude at the onset of reheating must take into account the effects 
of spacetime dynamics. This should be carried out self-consistently using the coupled semi- 
classical Einstein equation and matter-held equations, so that no ad hoc assumptions need 
be made about the effective equation of state and/or the relevant time scales involved. 

A full two-loop treatment of the unbroken symmetry mean-held dynamics of the O (N) 
held theory [which involves solving the nonlocal, integro-differential equations (5.24) and 
(5.25) of Ref. |IJ] includes multiparticle scattering processes, which provide a mechanism 
for reheating; but they are of a qualitatively different nature than the parametric resonance 
energy-transfer mechanism studied here. In addition, the nonlocal nature of the gap equation 
in the full two-loop analysis makes numerical solution of the coupled Einstein and matter 
equations difficult. In our model, multiparticle scattering occurs on a time scale T 2 which 
is significantly longer than the time scale T\ for parametric amplihcation of quantum huc- 


8 In recent work on galaxy formation from quantum fluctuations, Calzetta, Hu, and Matacz [55,56] 
report that A can be as high as ~ 10 -5 . 
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tuations. Therefore, in this model the leading-order, large -N (collisionless) approximation 
is sufficient for a study of parametric amplification of quantum fluctuations. In addition, 
realistic models of inflation invariably involve couplings of the inflaton to other fields, which 
provides additional mechanisms of energy transfer away from the inflaton mean field, and 
into its (or other fields’) quantum modes. 

The issues involved in a systematic study of the thermalization stage of post-inflationary 
physics are more complex. A quantum kinetic field theory treatment taking into account 
multiparticle scattering is required. The two-loop, 2PI effective action is the simplest and 
most generally applicable rigorous formalism which contains the leading-order multiparti¬ 
cle scattering processes. The leading-order, 1/N approximation is a collisionless subcase 
of the two-loop, 2PI effective action; it is employed in this study in order to obtain local 
dynamical equations which can be solved numerically, and is adequate for a study of para¬ 
metric amplification of quantum fluctuations. In addition, the growth of entropy must be 
understood within the context of a physically meaningful coarse graining of the full time- 
reversal-invariant quantum dynamics of the field theory. This, together with a first-principles 
analysis of the thermalization stage, is currently under investigation |fA| |87| . 
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TABLES 


Figures 

Hvo) 

A 

M P 

I\ 

o 

7 

IB 


2 x 10 7 

1 x 1(T 14 

1 x 10 14 

50.0 

1.7275 x 10 

PH 

1 


2 x 10 7 

1 x 1(T 14 

1 x 10 12 

50.0 

1.7275 x 10 

PP 

14 


2 x 10 7 

1 x 10~ 14 

6 x 10 10 

50.0 

1.0364 x 10 

15 

-19 

2 x 10 7 

1 x 1(T 14 

6 x 10 9 

50.0 

103.65 


TABLE I. Values of parameters for numerical solutions of Eqs. Q , ( El ), ( |3.30| ) in units 
where m = 1, and the corresponding figures in which the solutions are plotted. The -ff _1 (?7o) 
column is the initial inverse Hubble constant, which gives the initial time scale for cosmic expansion. 
Figs. 0H plot the resulting solutions. 


43 





























44 







































































































































































































45 





























































46 



























































































































































47 






































































































































































































































































































FIG. 14. Plot of p Q vs 77 , with M P /m = 6.0 x 10 10 . 
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